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O ! ABSTRACT 



^1- 

o 



We have carried out a high spectral resolution (A/AA ~ 6800 — 9700) line survey 
towards the Orion Klcinmann-Low (KL) cluster from 44-188 /im. The observations 
were taken with the Long Wavelength Spectrometer (LWS) in Fabry-Perot mode, on 
' board the Infrared Space Observatory (ISO). A total of 152 lines are clearly detected 

, and a further 34 features are present as possible detections. The spectrum is dominated 

■ by the molecular species H2O, OH and CO, along with [O 1] and [C 11] lines from PDR 

' or shocked gas and [O in], [N ill] lines from the foreground M42 H 11 region. Several 

isotopic species, as well as NH3, are also detected. HDO and ti30 + are tentatively 
Q ■ detected for the first time in the far-infrared range towards Orion-KL. A basic analysis 

l-j I of the line observations is carried out, by comparing with previous measurements and 

CZ2 . published models and deriving rotational temperatures and column densities in the 

case of the molecular species. 

Analysis of the [O 1] and [C 11] fine structure lines indicates that although a shock 
model can reproduce the observed [O 1] surface brightness levels, it falls short of 
the observed [C 11] level by more than a factor of 30. A PDR model can reproduce 
^ , the [O 1] 63.2 /im and [C 11] surface brightness levels within 35%, although over- 

predicting the LWS [O 1] 145.5 fxm emission by a factor of 2.7. The 70 water lines 
and 22 OH lines detected by the survey appear with mainly P Cygni profiles at the 
shortest survey wavelengths and with mainly pure emission profiles at the longest 
survey wavelengths. The emission and absorption velocity peaks of the water and 
OH lines indicate that they are associated with gas expanding in the outflow from 
the KL cluster. The estimated column densities are (2-5) xlO 14 cm~ 2 for H2O and 
(2.5-5.1) x 10 16 cm" 2 for OH. The 26 detected CO lines confirm the presence of three 
distinct components, with temperature and column density combinations ranging from 
660 K, 6 x 10 17 cm -2 to 360 K, 2 x 10 19 cm -2 . Comparison of the surface brightnesses 
and integrated fluxes of the CO lines measured in the 80-arcsec LWS beam with those 
measured previously by the Kuiper Airborne Observatory (KAO) in a 44 arcsec beam 
shows similar surface brightnesses in the different beams for the lowest-J CO lines and 
similar integrated fluxes in the different beams for the highcst-J CO lines, indicating 
that emission from the former lines has a uniform surface brightness within the LWS 
beam, while the latter lines originate from a region less than 44 arcsec in diameter. 
The complexity of the region requires more sophisticated models for the interpretation 
of all the line observations. 

Key words: infrared: ISM - ISM: molecules - ISM: individual (Orion) - surveys - 
line: identification - ISM: lines and bands 
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1 INTRODUCTION 

The formation of a star is often associated with outflows 
which can compress and heat the surrounding quiescent 
gas, altering the chemistry of the parent molecular cloud in 
which the star is forming. Many of the molecules responsible 
for cooling the gas in these molecular clouds, such as CO 
(J ^ 14), OH and H2O, emit at far-infrared wavelengths 
(> 40/xm) which are difficult to detect from the ground due 
to the constraints of the Earth's atmosphere. The launch 
of ESA's Infrared Space Observatory (ISO) in 1995 marked 
a revolution in our ability to observe these transitions, 
and in particular, opened a window for large-scale high 
resolution far-infrared spectroscopy of star forming regions 
such as Orion. In this paper we report on the first, and to 
date only, unbiased high spectral resolution spectral survey 
in the far-infrared waveband of the Kleinmann-Low (KL) 
region of Orion and discuss the implications of the results 
for the interpretation of the chemical evolution of the region. 

The Orion star-forming complex, at a distance of 450 
pc, is our nearest region of high mass star formation. It 
mainly consists of two giant molecular clouds: Orion molec- 
ular clouds 1 and 2 (OMC1 and OMC2), which were revealed 
via large scale mapping in the J=1^0 12 CO rotational tran- 
sition (Tucker, Kutner & Thaddeus, 1973). OMC1 is asso- 
ciated with the dark clouds L1640, L1641, L1647 (Lynds 
1962), extending over 6° southwards from the Orion nebula. 
The OMC1 ridge is associated with regions of prominent 
molecular emission, with the region near the Orion Neb- 
ula being the best studied example. OMC1 contains the 
visible nebula M42, a blister of hot, photo-ionised lumi- 
nous gas around the hot Trapezium stars (see Figure^ and 
also contains a number of IR-emitting regions such as the 
Kleinmann-Low nebula, which is the brightest far-infrared 
region in the complex (Kleinmann & Low 1967). The KL 
nebula includes an infrared cluster of massive stars at an 
early evolutionary stage, e.g. the Becklin-Neugebauer ob- 
ject (BN, Becklin & Neugebauer 1967) and IRc2, which 
dominate the mid-infrared radiation (Downes et al. 1981). 
There are at least two molecular outflows, one associated 
with material ejected at hundreds of km s _1 forming a bipo- 
lar cone of molecular fingers, whose axis is perpendicular to 
the NH3 emission associated with the hot core (Wilson et 
al. 2000), and a second low- velocity outflow first detected 
by H2O maser observations (Genzel et al. 1981). IRc2 was 
originally believed to be the major source of the KL lumi- 
nosity (L~ 10 s L Q , Genzel & Stutzki, 1989). However, the 
infrared observations of Dougados et al.(1993) resolved IRc2 
into four components which may not even be self-luminous, 
and radio studies located the origin of the outflow as being 
offset by 0.5 arcsec from IRc2-A; at the position of radio 
source I (Menten & Reid 1995, Gezari et al. 1998, Chandler 
& Greenhill 2002). 

Within the KL region, molecular emission arises from sev- 
eral physically distinct regions: 

• The Hot Core is composed of very dense material with 
n ~ 10 7 cm -3 and T ~ 200K in clumps located ~ 2" south 
of IRc2, 1" offset from radio source I (Wright et al. 1996, 
Schilke et al. 2001). 

• The Compact Ridge is a compact region (14") of dense 
gas (n ~ 10 6 cm -3 ) surrounded by outflows. 



• The Plateau; a region of inhomogeneous density con- 
taining both low-velocity (~ 18 km s _1 ) and high-velocity 
(~ 100 km s _1 ) flows with an approximate diameter of 50" 
(Blake et al. 1987). 

• Finally the whole region is embedded in the Extended 
Ridge, a region of quiescent and cooler gas from which 
many narrow lines (Av ~ 3-5 km s _1 FWHM) are detected 
(Greaves & White 1991). 

The proximity and complexity of the KL region have made it 
the target of many astronomical studies, and it is one of the 
best studied examples of the interaction between massive 
stars in their earliest stages and their parental molecular 
cloud. In this paper we report the results of an unbiased 
spectral survey between 44-188/im (1612 to 6813 GHz) at 
a resolving power of ~ 8000, which represents the first sys- 
tematic study of the molecular spectrum in this wavelength 
range. A number of the detected lines from this spectral sur- 
vey and from a series of targeted observations have already 
been reported (e.g. the OH lines; Goicoechea et al. 2006, 
some H2O lines; Harwit et al. 1998, Cernicharo et al. 1999, 
and some CO lines Sempere et al. 2000), however, this is the 
first analysis of the entire spectrum. 

Our spectral survey provides an insight into the physical 
conditions and dynamical processes in the Orion KL region 
and complements studies performed at lower frequencies, 
e.g. from 72 to 91 GHz (Johansson et al. 1984), 70 to 115 
GHz (Turner 1989), 150 to 160 GHz (Ziurys & McGonagle 
1993), 215 to 247 GHz (Sutton et al. 1985), 247 to 263 GHz 
(Blake et al. 1986) 325 to 360 GHz (Schilke et al. 1997), 455 
to 507 GHz (White et al. 2003), 607 to 725 GHz (Schilke et 
al. 2001) and 795 to 903 GHz (Comito et al. 2005). 
The aim of this work is to present the data, describe its re- 
duction and calibration and discuss some basic results. Sec- 
tion 2 describes the observations and data reduction, includ- 
ing a brief summary of the instrumental operation. Section 
3 presents the results, while the individual species are dis- 
cussed in Section 4. Our results are summarised in Section 
5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Instrumental operation 

The Long Wavelength Spectrometer (LWS; Clegg et al. 
1996) was one of the four instruments on board the 
Infrared Space Observatory (ISO), which operated between 
November 1995 and April 1998, when the superfluid helium 
used to cool the instruments ran out. The LWS covered 
the spectral range between 43 and 197^im at medium 
resolution (A/AA ~ 150-200) using a diffraction grating 
and at high resolution (A/AA ~ 6800-9700) with one of 
two Fabry-Perot (FP) etalons, order sorted by the grating. 
Four different observing modes were available to users, 
via the so-called 'Astronomical Observation Templates' 
(AOTs), which allowed the observer a choice of wavelength 
ranges, sampling intervals and exposure times. These 
observing modes were: medium resolution, full wavelength 
range spectrum (AOT L01), medium resolution and narrow 
band photometry, selected wavelength range spectrum 
(AOT L02), high-resolution, full wavelength range (AOT 
L03) and high-resolution, narrow wavelength range (AOT 
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Figure 1. Colour composite mosaic image of the central part 
of the Orion Nebula, M42 , based on 81 images obtained with 
the infrared multi-mode ISAAC instrument on the ESO Very 
Large Telescope (VLT) at the Paranal Observatory. The Trapez- 
ium stars are seen near the centre (ESO PR Photo 03a/01 2001). 



L04). Also, a parallel mode was carried out while an AOT 
was active with another instrument and a serendipity mode 
was carried out during the slews with no prime instrument. 
Since a full spectral scan with the LWS FP (AOT L03) 
took a long time, a significant part of the spectrum was 
recorded at high spectral resolution for only four objects: 
Orion BN/KL, Sgr B2 (Polehampton 2002, Polehampton 
et al. 2006 in prep.), Sgr A and Jupiter. 

In L03 mode, data for a given wavelength range were 
recorded by scanning the FP gap over a small range with the 
grating held stationary, thus shifting the wavelength of the 
selected FP order across the grating response function. To 
record data for another wavelength range the grating was 
moved by a portion of the grating resolution element and 
the FP scanned again: the data set recorded at each grating 
setting is termed a mini-scan. 

During each observation the LWS FP and grating settings 
were optimised for the detector whose band pass filter in- 
cluded the wavelength range of interest. This was denoted 
as the 'prime' detector. However, all ten LWS detectors 
recorded data simultaneously in their own spectral ranges; 
the other detectors are termed 'non-prime'. 

2.2 Observations 

The observations were carried out between September 
1997 and April 1998. The dataset consists of 26 individual 
observations making up a total of 27.9 hours of ISO LWS 
observing time in L03 mode, 16 observations making up 
13.1 hours in L04 mode and 1 observation in the lower 
resolution L01 grating mode. The instrumental field of view 
for all L03 observations was either centred on a position 
offset by 10.5" from the BN object (which is at 05 h 35 m 



14.12 s - 05° 22' 22.9" J2000), or a position offset by 5.4" 
from IRc2 (which is at 05 h 35 m 14.45 s - 05° 22' 30.0" 
J2000) while most of the L04 observations were centred 
on IRc2 (see Table □ and Figure ^J. The LWS beam had 
a diameter ^ 80" (Gry et al. 2003). Table [T] lists the 
observations: column one is the Target Dedicated Time 
identification number (TDT). The TDT is an eight digit 
number unique to each observation; the first three digits 
are the ISO revolution number, the next three digits are 
the sequence number of the observation in the revolution 
and the last two numbers were assigned by the observer. 
Column two lists the wavelength range covered by the prime 
detector. The length of the observation in seconds can be 
found in column three while the observation date (d.m.y) 
is in column four. The J2000 coordinates of the pointing 
position for each observation can be found in columns five 
and six respectively. 



2.3 Data reduction 

Processing of the LWS FP data was carried out using the 
Offline Processing (OLP) pipeline and the LWS Interactive 
Analysis (LIA) package version 10. The basic calibration 
is fully described in the LWS handbook (Gry et al. 2003). 
Further processing, including dark current optimisation, de- 
glitching and removal of the LWS grating profile was then 
carried out interactively using the LWS Interactive Analysis 
package version 10 (LIA10; Lim et al. 2002) and the ISO 
Spectral Analysis Package (ISAP; Sturm 1998). The reduc- 
tion process was performed in the same way as that for the 
ISO L03 Sagittarius B2 data (Polehampton et al. 2006 in 
prep) with the only exception being the mini-scan shape 
removal (see below). The dataset contains data from both 
prime and non-prime detectors. 

Detector dark currents were determined during the mission 
by three different methods (Swinyard et al. 1998). Table |2] 
lists the nominal dark current, which is the average of these 
measurements. However, scattering of light into the detec- 
tors may have occurred during observations of extended or 
off-axis sources. This could contribute an additional compo- 
nent to the dark signal which varied from source to source 
and is not easy to distinguish from the detector dark current 
or the source signal. A method to determine both the dark 
current and stray light corrections for a particular observa- 
tion was developed by Polehampton (2002) using non-prime 
data from the observation of Sagittarius B2 (see also Pole- 
hampton et al. 2006). The method is based on the fact that 
for some mini-scans the transmitted orders from the FP do 
not correspond to the wavelength range transmitted by the 
grating. In these circumstances no light was transmitted and 
the detector will have recorded only the dark current and the 
stray light that finds its way to the detectors from outside 
of the main beam - the dark signal. 

This method was used on the Orion data and dark signal 
corrections were applied via the graphical interface available 
in the LIA interactive reduction package. Table |2] lists the 
values of the dark signal that were subtracted from the ob- 
servations for each detector compared to the 'nominal' dark 
current. 

Basic removal of the mini-scan shape was also performed 
interactively in the LIA. Any remaining mini-scan residual 
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TDT Wavelength(/im) 


lcngth(s) 


Date 


J2000 RA 


J2000 Dec 


AOT = L03 


66302402 


47- 


-52 


5586 


9.9.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


5 


66302406 


47 


52 


5586 


9.9.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


5 


66002003 


52 


57 


4434 


6.9.97 


5> 


35 m 


14 s 


.17 


-5° 


22' 


33" 


4 


66002007 


52 


57 


4434 


6.9.97 


5 1 


35 m 


11 s 


.17 


-5° 


22' 


33" 


1 


70101704 


57 


63 


5726 


17.10.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


70101708 


57 


63 


5726 


17.10.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


70001205 


63 


70 


5641 


16.10.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


70001209 


63 


70 


5642 


16.10.97 


5> 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


69901510 


70 


73 


2592 


15.10.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


69901514 


70 


73 


2591 


15.10.97 


5 h 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


69901312 


77 


-81 


4106 


15.10.97 


5 1 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


69901316 


77 


-81 


4106 


15.10.97 


5' 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


69901413 


81 


85 


3842 


15.10.97 


rh 

□ 


n r m 
OO 


14 s 


.17 


r o 
-O 


22' 


on" 
66 


8 


69901417 


81 


85 


3842 


15.10.97 




35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


70001105 


85- 


89 


2530 


15.10.97 


5' 


35 m 


14 s 


.16 


-5° 


22' 


33" 


7 


70001127 


85 


89 


2530 


16.10.97 


5 1 


35 m 


14 s 


.17 


-5° 


22' 


33" 


8 


87301008 


99 


104 


2924 


6.4.98 


5 1 


35 m 


14 s 


.15 


-5° 


22' 


33" 


1 


87301030 


99 


104 


2924 


6.4.98 


5 h 


35 m 


If 


.15 


-5° 


22' 


33" 


4 


87300909 


104- 


-109 


2748 


6.4.98 


5 h 


35 m 


14 s 


.15 


-5° 


22' 


33" 


1 


87300931 


104- 


-109 


2748 


6.4.98 


5' 


35 m 


11 s 


.15 


-5 C 


22' 


33" 


■1 


82702210 


109- 


-115 


3226 


19.2.98 


5 h 


35 m 


ir 


.15 


-5° 


22' 


33" 


5 


82702232 


109- 


-115 


3226 


20.2.98 


5 h 


35 m 


14 s 


.15 


-5° 


22' 


33" 


5 


84101911 


115- 


-121 


3140 


5.3.98 


5' 


35 m 


14 s 


.15 


-5° 


22' 


33" 


■1 


84101933 


115 


-121 


3140 


5.3.98 


5 1 


35 m 


14 s 


.15 


-5° 


22' 


33" 


5 


82602316 


147- 


-154 


3830 


18.2.98 


5 1 


35 m 


14 s 


.15 


-5° 


22' 


33" 


5 


82602338 


147- 


-154 


3830 


18.2.98 


5 h 


35 m 


If 


.15 


-5° 


22' 


33" 


5 


AOT = L04 


83101203 


51- 


101 


2166 


23.02.98 


5 h 


35 m 


14 s 


.45 


-5° 


22' 


30" 





69000714 


51 


125 


6918 


06.10.95 


5 1 


35 m 


14 s 


.21 


-5° 


22' 


23" 


5 


82901207 


51 


75 


1818 


22.02.98 


5 h 


35 m 


14 s 


.45 


-5° 


22' 


30" 





82901101 


68 


129 


3858 


21.02.98 


5 h 


35 m 


14 s 


.45 


-5° 


22' 


29" 


5 


69602606 


70 


124 


2334 


12.10.97 


5' 


35 m 


14 s 


.47 


-5° 


22' 


29" 


8 


69602605 


73 


180 


2028 


12.10.9 


5 1 


35 m 


14 s 


.47 


-5° 


22' 


29" 


8 


82901206 


75 


91 


1772 


21.02.98 


5 h 


35 m 


14 s 


.45 


-5° 


22' 


30" 


1 


83101202 


83 


180 


2282 


23.02.98 


5 h 


35 m 


14 s 


.45 


-5° 


22' 


30" 


1 


70101216 


84- 


179 


3028 


17.10.97 


5 h 35 m 


14 s 


47 


-5° 


22' 


30" 


4 


83101201 


96 


180 


2416 


23.02.98 


5 1 


35 m 


14 s 


.45 


-5° 


22' 


30" 


1 


69602317 


100- 


-180 


1980 


12.10.97 


5 h 35 m 


14 s 


47 


-5° 


22' 


30" 


1 


82901204 


105- 


-146 


1630 


21.02.98 


5* 


35 m 


14 s 


.45 


-5° 


22' 


30" 


1 


83101301 


119- 


-120 


3380 


23.02.98 


5> 


35 m 


14 s 


.49 


-5° 


22' 


30" 


9 


68701716 


119 


-165 


7661 


03.10.97 


5 1 


35 m 


14 s 


.21 


-5° 


22' 


23" 


5 


69602318 


119 


-179 


2492 


12.10.97 


5 h 


35 m 


14 s 


.48 


-5° 


22' 


30" 


1 


68701301 


148- 


-152 


1294 


03.10.97 


5 h 


35 m 


14 s 


.47 


-5° 


22' 


29" 


8 



Table 1. List of all L03 and L04 observations centred on Orion IRc2-BN/KL analysed in this work. Observations are sorted by wavelength 
range covered by the prime detector. Note that for L04 spectra, only certain targeted observations were covered in the quoted wavelength 
range. 



shape is due to incorrect substraction of the grating response 
profile as a consequence of a positional uncertainty in the ab- 
solute grating position. The effect can be seen as periodic 
structures in regions where mini-scans overlap. These fea- 
tures are very difficult to remove being still present in the 
final LI A product. Therefore, a further step was performed 
for the Orion data by using smooth functions and removing 
the residuals individually for every mini-scan (see also fol- 
lowing section). 

The Orion KL region was also observed using the grating 
in L01 mode in revolution 695. Due to the strength of the 
source, the grating data had to be corrected for the effects of 



detector saturation (Leeks et al. 1999). Channel fringing and 
spurious features introduced by the Relative Spectral Re- 
sponse Function (RSRF) were also corrected by using the 
Highly Processed Data Products (HPDP) available in the 
ISO data archive (Lloyd et al. 2003). The final L01 product 
has a superior flux calibration accuracy compared to the FP 
data and was therefore used to derive the continuum level. 
As part of the normalisation process, FP data were first 
smoothed to the same resolution as the grating, using a bin 
size of 0.07 fj,m for the SW detectors and 0.15 pm for the 
LW detectors. Finally, each FP spectrum was corrected to 
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Detector Nominal FPS FPL 



dark current (x lCT 16 A) (x 1CT 16 A) (x lCT 16 A) 



SW1 


4.96 


± 


0.55 














SW2 


2.08 


± 


0.43 


6.10 


± 


0.20 








SW3 


2.20 


± 


0.21 


4.97 


± 


0.19 








SW4 


1.18 


± 


0.34 


2.41 


± 


0.11 


5.48 


± 


0.49 


SW5 


1.56 


± 


0.24 


2.11 


± 


0.06 


5.75 


± 


0.56 


LW1 


2.50 


± 


0.29 


3.80 


± 


0.15 


4.20 


± 


0.07 


LW2 


0.07 


± 


0.27 


1.03 


± 


0.05 


4.64 


± 


0.06 


LW3 


0.53 


± 


0.39 


2.01 


± 


0.42 


4.83 


± 


0.35 


LW4 


1.76 


± 


0.42 


4.41 


± 


0.65 


5.28 


± 


0.28 


LW5 


1.21 


± 


0.25 


1.49 


± 


0.14 


1.47 


± 


0.06 



Table 2. Nominal dark currents adopted in the standard pipeline processing for OLP version 10 (column one) and averaged values of 
the dark current and stray light correction determined for the Orion dataset for the two FPs. 



the Local Standard of Rest velocity (LSR) using the Starlink 
routine RV (Wallace & Clayton 1996). 

2.4 Spectrum reliability, signal to noise achieved 

After the dark current optimization and mini-scan shifting 
many spurious features were observed in the data. These 
structures can be easily recognized and divided into two 
different types: features due to mini-scans shifting, and 
features transferred via the Relative Spectral Response 
Function (RSRF) in the calibration stage. The first are 
introduced by the fact that each observation is built up 
of mini-scans and interactively adjusted as part of the 
reduction process. 

This residual shape is difficult to remove, especially for 
short wavelength detectors and an alternative method was 
used to minimize transferred features. For each mini-scan 
a smooth function was fitted and used to remove the 
mini-scan shape. This method significantly improved the 
signal to noise achieved but introduced an error due to 
the fact that each observation is divided by the flux level 
achieved in each particular mini-scan. This means that the 
same lines recorded in different observations could present 
a slight difference in flux level in the final reduced file. 
However, the difference in flux when the same spectral 
range was observed by two adjacent detectors was found to 
be relatively small (see Figure [3J , introducing errors of less 
than 20%. 

The second group of features are transferred from the 
calibration via the RSRF and are probably caused by 
transient effects in the detectors. They can be seen as large 
scale features. 

The overall signal to noise achieved in the survey as a 
function of wavelength is shown in Figure 2] The inclusion 
of non-prime data fills in some gaps in spectral coverage 
that were present in the prime data as well as overlapping 
in some places with the prime observations. This improves 
the signal to noise achieved by a factor that depends on the 
detector. Polehampton et al. (2006 in prep.) give a detailed 
explanation of the effect of including non-prime data on the 
overall signal to noise, based on the throughput of the LWS 



FP etalons. 

High spectral resolution can provide detailed information as 
a function of velocity for atomic and molecular transitions 
as well as unique information about the temperature and 
density regimes in which the transitions occur. Information 
on dynamical processes can be inferred by comparing line 
shapes; this can be easily seen for Orion KL, with H2O and 
OH lines showing pure absorption, P-Cygni and pure emis- 
sion profiles (see Section 4). 

The advantages of high resolution spectroscopy can be seen 
by comparing the same wavelength region observed in the 
low resolution (L01) and high resolution (L03) mode. Figure 
shows a section of the low spectral resolution observation 
(L01 mode) compared with the high spectral resolution (L03 
mode) scan of the same region. At least 80% of the lines de- 
tected in the L03 scan are missed in the lower resolution 
spectrum. 

2.5 The LWS Beam profile 

The LWS beam profile is the convolution of the telescope 
point spread function (PSF) with the aperture of each de- 
tector. Ideally, the telescope PSF should be an Airy profile 
but the central obscuration, secondary supports and any op- 
tical imperfections redistribute power from the core of the 
profile to the Airy rings (Gry et al. 2003). The latest model 
of the telescope PSF (Gry et al. 2003) includes the effects 
of the central obscuration and its supporting structure, and 
indicates that the power in the Airy rings is increased and 
that the wings of the profile contain 2D asymmetric struc- 
ture. The asymmetry introduced into the profile is due to 
the three-legged secondary support. The aperture was as- 
sumed to be circular with a top hat profile (Gry et al. 2003, 
Lloyd et al. 2003). Table [3] lists the effective diameter of the 
aperture of each detector used to derive the effective solid 
angle. The solid angle is defined as 27r(l — cos 9), where 8 is 
the effective beam angular radius. 

2.6 Line flux measurements 

The shape of the FP spectral response function was deter- 
mined from dedicated calibration observations of narrow, 
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Figure 2. Approximate LWS observation center (cross points) 
at the positions of IRc2 and 10.5" offset from the BN object on 
a false colour-image from Lis et al. (1998) showing the 350/^m 
continuum emission observed with SHARC. (Courtesy of D. Lis). 
The circles show the LWS beam positions. 



unresolved spectral lines in the spectra of planetary nebula. 
These observations show that the response function shape 
can be accurately described by the predicted (Hecht and Za- 
jac, 1974) Airy profile. The observed line profile is then the 
convolution of the source intrinsic profile with the Airy func- 
tion. Line measurements were done by fitting a Lorentzian 
profile using the interactive line-fitting routine in ISAP. The 
Lorentzian function was found to approximate the line shape 




50 100 150 200 

Microns 



Figure 3. All L03 prime observations (upper plot) and non-prime 
observations (down plot) together with a continuum fit to the 
grating observations. Note that for wavelengths ^120 (im the 
spectrum is mainly covered by non-prime data. 



150 ■- 




1 SO 



Wavelength (/i-m) 

Figure 4. The signal to noise achieved for the L03 prime-detector 
spectra. Peaks in the curve represent the spectral ranges with 
the highest data reliability. The vertical dotted line divides the 
spectral ranges observed by the short wavelength FP (FPS) from 
those observed by the long wavelength FP (FPL) 



very well, however, for unidentified features and lines with 
poor signal, a Gaussian fitting function was selected. Data 
were averaged over detectors and scans. The ISAP fitting 
procedure gives the line flux, line centre, FWHM, and an 
estimate of the statistical error on the line flux. 
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Figure 5. Comparison of a section of the low resolution L01 
spectrum of Orion KL (top plot) with the high resolution L03 
spectrum (middle plot). The bottom plot shows a direct compar- 
ison of the two spectra in a narrower spectral range. 



Detector 


Wavelength 
range (/tm) 


Effective 
radius (") 


Effective 
solid angle (sr) 


SW1 


43-50.5 


39.4 


1.15 X 


io- 7 


SW2 


49.5-64 


42.3 


1.32 x 


io- 7 


SW3 


57-70 


43.5 


1.40 x 


IO" 7 


SW4 


67-82 


40.9 


1.23 x 


IO" 7 


SW5 


76-93 


39.5 


1.15 x 


10~ 7 


LW1 


84-110 


38.6 


1.10 x 


io- 7 


LW2 


103-128 


38.9 


1.12 x 


IO" 7 


LW3 


123-152 


35.5 


9.30 x 


IO" 8 


LW4 


142-171 


34.7 


8.89 x 


10" 8 


LW5 


167-197 


33.2 


8.14 x 


IO" 8 



are also detected. Table [I] and Figure |S| show the identified 
lines. 



4 DISCUSSION OF INDIVIDUAL SPECIES 

A basic analysis and discussion of the individual detected 
species is presented in the following sub-sections. 
From the analysis of the atomic forbidden lines, we derive 
physical parameters such as temperatures and densities by 
comparing our measurements with models. When conditions 
are appropriate, we use the 'rotation diagram' method to de- 
termine molecular parameters such as rotation temperatures 
and column densities. A problem with using the rotational 
diagram method is that it underestimates the total column 
density if the energy distribution is far from Local Thermo- 
dynamic Equilibrium (LTE). In order to decide whether the 
use of this method is appropriate, we estimated line opti- 
cal depths using the statistical equilibrium radiative trans- 
fer code RADEX. (Schoier et al. 2005). RADEX is a one- 
dimensional non-LTE radiative transfer code that uses the 
escape probability formulation assuming an isothermal and 
homogeneous medium without large-scale velocity fields. 



4.1 Atomic forbidden Lines 

Table |S| lists several detected forbidden transitions observed 
by the survey: [O ill] 52 fim, 88 fim, [N ill] 57 /im, [N n] 
122 fim, [O i] 63 fim, 145 fim and [C n] 158 ^m. The flux 
errors were deduced from Lorentzian profile fitting (except 
the [N n] 122 fxm line) and do not include systematic errors 
caused by uncertainties in the flux calibration, see Section 
2.4. The [N n] 122 /im line was fitted using a Gaussian profile 
due to its width and low intensity. This broadness is due to 
hyperfine components in the transition as can be seen in the 
detection of the [N n] 122 /im line in Sgr B2 (Polehampton 
et al. 2006, in prep.). 



4-1.1 Lines from ionised gas: [O ill], [N ill] and [N n] 

The [O ill] 52, 88 (im, [N III] 57 fim and [N n] 122 /im lines 
originate from the foreground M42 H II region which is ex- 
cited by the 9 1 Ori OB stars, known as 'The Trapezium' (see 
Figure 0; these lines provide an important tool for probing 
physical conditions and elemental abundances in the ionised 
gas (Simpson et al. 1986, Rubin et al. 1994). 



Table 3. The effective beam radius and effective solid angle sub- 
tended by each detector. Differences in effective apertures are due 
to asymmetries introduced into the profile by the three-legged 
secondary instrumental support (Gry et al. 2003) 



3 DATA ANALYSIS AND RESULTS 

Line identification was carried out based on the JPL line 
catalogue (Pickett et al. 1998) and the molecular catalogue 
created and maintained by one of us (JC), which includes 
more than 1200 molecular species. The spectrum is domi- 
nated by molecular lines of CO, OH and H2O and by for- 
bidden lines of [O ill], [O 1], [N ill], [N 11] and [C 11]. Isotopic 
variants and the lowest transitions of HDO, NH3 and HaO + 



4-1.2 Electron densities from the [O ill] lines 

The intensity ratio of the [O ill] 52 /im and 88 /im lines 
provides an electron density diagnostic. The line ratio has 
almost no dependence on electron temperature because of 
the very low excitation energies of the transitions. The mea- 
sured line fluxes give an 88.4 /im/51.8 /im flux ratio of ~ 
0.41 ± 0.01. Adopting the same atomic parameters as Liu et 
al. (2001), we obtain an electron density iV e (0 m)=830±i2o 
cm -3 calculated for T e =8000 K. 

Another important electron density indicator is the [S ill] 
33.5/18.7 p,m ratio. 
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Table 4: Line measurements obtained from L03 and L04 high resolution observa- 
tions towards the Orion KL region. Line fitting measurements were performed 
using a Lorentzian function. Quoted errors are based on goodness of fit esti- 
mations. Note that ( x ) corresponds to a tentative detection with flux below 
3<t; ( 2 ) are blended lines and ( 3 ' are lines detected in the lower resolution L01 
spectrum where no velocity peak information is given. 



Species Transition Rest Flux Flux Peak Peak 

Wavelength absorption emission absorption emission 

(/jm) (10~ 1T W cm -2 ) (10~ 17 W cm -2 ) (km s _1 ) (km s _1 ) 



P-H 2 


5 42 - 431 


44.195 


0.48 ± 0.091 


0.39 ± 0.20 


-28.9 ± 5.5 


53.6 ± 27.4 


o-H 2 


523 - 4i4 


45.111 


1.11 ± 0.016 


0.39 ± 0.091 


-22.2 ± 2.1 


39.7 ± 9.2 


P-H 2 


3.31 - 2o2 


46.484 


1.05 ± 0.12 




-10.6 ± 2.1 




P-H 2 


IO19 - 928 


47.039 


0.19 ± 0.041 


0.49 ± 0.19 


-3.1 ± 1.1 


23.6 ± 9.1 


o-H 2 


532 - 4 23 


47.972 


1.13 ± 0.14 


0.55 ± 0.21 


-16.3 ± 2.1 


38.9 ± 14.8 


OH 


2 n 1/2 - 2 n 3/2 .1=5/2- - 5/2+ 


48.704 


1.09 ± 0.13 




-12.4 ± 1.4 




OH 


2 n 1/2 - 2 n 3/2 j=5/2+ - 5/2- 


48.817 


1.34 ± 0.14 




-22.8 ± 2.8 




P-H 2 


440 — 331 


49.281 


0.71 ± 0.13 




-9.2 ± 1.6 




o-H 2 


441 — 330 


49.336 


1.82 ± 0.26 


0.63 ± 0.14 


-17.1 ± 2.4 


48.1 ± 10.6 


o-H 2 


634 - 5 23 


49.390 




0.99 ± 0.21 




44.1 ± 9.3 


P-H 2 


928 - 8 17 


50.634 




0.19 ± 0.02 




35.4 ± 3.7 


CO' 1 ' 


52 - 51 


50.887 




0.71 ± 0.23 




1.1 ± 0.4 


[OIII] 


3 P 2 - 3 Pi 


51.815 




37.1 ± 0.50 




7.2 ± 0.1 


P-H 2 


533 — 4 2 2 


53.137 


0.24 ± 0.04 


0.41 ± 0.03 


-32.1 ± 5.4 


36.5 ± 2.7 


OH 


2 n 1/2 - 2 n 3/2 j= 3/2+ - 3/2- 


53.261 


2.55 ± 0.32 


1.10 ± 0.62 


-20.7 ± 2.5 


14.5 ± 8.1 


OH 


2 n 1/2 - 2 n 3/2 j= 3/2- - 3/2+ 


53.351 


2.45 ± 0.26 


0.26 ± 0.13 


-14.2 ± 1.5 


45.1 ± 19.9 


P-H20 


431 - 322 


56.324 


1.09 ± 0.11 


0.59 ± 0.25 


-5.1 ± 0.5 


44.9 ± 18.9 


[NIII] 


P3 /2 ~ ^~*l/2 


57.329 




5.97 ± 0.12 




2.9 ±0.1 


P-H20 


4 22 — 3i3 


57.636 


2.15° ± 0.24 


0.52 ± 0.19 


-15.5 ± 1.7 


48.9 ± 17.8 


CO 


45 - 44 


58.547 




0.12 ± 0.037 




6.7 ± 2.1 


o-H 2 


432 — 3 2 i 


58.698 


1.57 ± 0.13 


1.80 ± 0.14 


-29.9 ± 2.5 


37.3 ± 2.9 


HDO 


330 — 2n 


59.928 




0.15 ± 0.013 




37.5 ± 3.2 


P-H 2 


431 - 4 4 


61.808 




0.52 ± 0.02 




43.5 ± 1.7 


HDO 


lOes - IO56 


62.231 




0.02 ± 0.006 




4.5 ± 1.3 


[OI] 


3 Pi - 3 P 2 


63.184 




65.33 ± 0.54 




14.6 ± 0.1 


P-H 2 


818 — ^07 


63.322 




0.30 ± 0.12 




30.3 ± 12.1 


o-H 2 


§08 - 7l7 


63.457 




0.27 ± 0.045 




41.5 ± 6.9 


OH 


2 n 3/2 j= 9/2- - 7/2+ 


65.132 


0.44 ± 0.091 


0.84 ± 0.12 


-22.9 ± 4.7 


29.6 ± 4.2 


o-H 2 


625 - 5l4 


65.166 




1.09 ± 0.13 




26.4 ± 3.1 


OH 


2 n 3/2 j= 9/2+ - 7/2- 


65.279 


0.39 ± 0.10 


0.59 ± 0.20 


-15.6 ± 3.9 


34.9 ± 11.8 


o-H 2 


7l6 ~~ 625 


66.092 




0.78 ± 0.16 




30.3 ± 6.2 


o-H 2 


330 ~~ 2 2 i 


66.437 


1.51 ± 0.18 


1.39 ± 0.28 


-22.2 ± 2.6 


40.5 ± 8.1 


P-H 2 


331 - 2 2 q 


67.089 


1.40 ± 0.21 


1.07 ± 0.20 


-12.7 ± 1.9 


42.5 ± 7.9 


H 2 8 ° 


330 ~~ 2 2 i 


67.191 


0.19 ± 0.031 




-12.9 ± 2.1 




0-H2O' 1 ' 


330 ~~ 3q3 


67.269 


0.33 ± 0.11 


0.64 ± 0.11 


-22.5 ± 7.5 


32.5 ± 5.6 


CO 


39 - 38 


67.336 




0.48 ± 0.14 




12.4 ± 3.6 


CO 


38 - 37 


69.074 




0.45 ± 0.14 




15.5 ± 4.8 


CO 


37 - 36 


70.907 




0.70 ± 0.061 




10.1 ± 0.9 


P-H 2 


5 2 4 - 4 13 


71.066 


0.15 ± 0.039 


0.89 ± 0.059 


-22.5 ± 5.8 


28.5 ± 1.8 


OH 


2 n 1/2 J= 7/2- - 5/2+ 


71.171 




0.33 ± 0.049 




28.2 ± 4.2 


OH 


2 n 1/2 J= 7/2+ - 5/2- 


71.218 




0.78 ± 0.12 




27.9 ± 4.3 


P-H 2 


7l7 - 606 


71.539 




1.39 ± 0.13 




22.3 ± 2.1 


o-H 2 


7o7 - 616 


71.946 




1.42 ± 0.28 




26.7 ± 5.2 


CO 


36 - 35 


72.843 




0.62 ± 0.13 




1.6 ± 0.3 


HDO' 1 ' 


9 6 3 - 1046 


74.792 




0.07 ± 0.02 




19.5 ± 5.6 


CO 


35 - 34 


74.890 




3.46 ± 0.40 




0.1 ± 0.01 


o-H 2 


3 2 i - 2 12 


75.380 


0.88 ± 0.25 


5.67 ± 0.21 


-31.5 ± 8.9 


28.5 ± 1.1 


o-H 2 


551 - 5 42 


75.779 




0.32 ± 0.04 




19.2 ± 2.4 


H 2 s O 


321 - 2 i 2 


75.868 


0.08 ± 0.01 




-10.1 ± 0.6 




o-H 2 


550 - 5 41 


75.908 




0.27 ± 0.01 




20.1 ± 0.7 


CO 


34 - 33 


77.059 




0.60 ± 0.09 




4.9 ± 0.7 


o-H 2 


4 2 3 - 3i 2 


78.741 




6.18 ± 0.17 




26.9 ± 0.7 


P-H 2 


615 - 5 2 4 


78.928 




0.43 ± 0.11 




16.5 ± 4.2 


OH 


2 n 3/2 j= 1/2- - 3/2+ 


79.118 




3.89 ± 0.16 




16.4 ± 0.7 


OH 


2 n 3/2 j= 1/2+ - 3/2- 


79.182 




3.73 ± 0.38 




19.1 ± 1.9 


CO 


33 - 32 


79.359 




2.32 ± 0.20 




3.7 ± 0.3 


CO 


32 - 31 


81.806 




1.52 ± 0.09 




10.8 ± 0.6 


o-H 2 


616 - 5o5 


82.030 




2.78 ± 0.13 




22.4 ± 1.1 


o-H 2 


836 - 8 2 7 


82.974 




0.27 ± 0.06 




15.1 ± 3.2 


P-H 2 


606 - 5i5 


83.283 




1.29 ± 0.08 




31.1 ± 1.8 


OH 


2 n 3/2 J= 7/2+ - 5/2- 


84.420 




1.72 ± 0.13 




30.8 ± 2.3 


O-H20' 1 ' 


1148 - IO55 


84.514 




0.15 ± 0.10 




2.1 ± 1.4 


OH 


2 n 3/2 J= 7/2- - 5/2+ 


84.597 


1.13 ± 0.38 


1.34 ± 0.26 


-21.4 ± 7.9 


36.2 ± 7.1 


o-H 2 (2) 


7 16 - 7 07 


84.766 




0.47 ± 0.05 




19.2 ± 1.9 


H^O* 1 ' 2 ' 


7 16 - 7 07 


84.778 










CO 


30 - 29 


87.190 




1.58 ± 0.11 




2.8 ± 0.2 


[OIII] 


3 Pl- 3 Po 


88.356 




13.2 ± 0.29 




9.8 ± 0.1 


o-H 2 


3 22 - 2 n 


89.988 




2.17 ± 0.17 




29.1 ± 2.3 


CO 


29 - 28 


90.163 




1.58 ± 0.14 




5.1 ± 0.4 



Tentative detections Blended lines L01 detection 
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Table 4 - Continued 



Species 



Transition 



Rest 
Wavelength 



Flux 
absorption 



Flux 

emission 



Peak 
absorption 



Peak 
emission 







f //ml 


(10- 17 W cm -2 ) 


(io- 17 w 


cm 2 ) 


(km s -1 ) 


(km 


s-i) 


p-H a O 


1157 - 1064 


90.643 




0.16 







03 




42.5 ± 8.0 


o-H 2 


643 - 634 


92.811 




0.49 


± 





00 




8.9 


± 1.1 


CO 


28 - 27 


93.349 




1.89 


± 





11 




10.8 


± 0.6 


o-H 2 


7 35 - 7 26 


93.379 




0.36 


± 





05 




20.5 


± 2.9 


P-H 2 


542 - 533 


94.206 




0.62 


± 





02 




14.6 


± 0.5 


o-H 2 


625 - 616 


94.643 




0.52 


± 





06 




22.7 


± 0.3 


o-H 2 


441 - 4 32 


94.703 




0.72 


± 





06 




27. 1 


± 2.2 


P-H 2 


5 15 - 4 4 


95.626 




1.42 


± 





10 




27.4 


± 1.3 


P-H 2 


4 4 1 - 4 32 


95.883 




0.67 


± 





06 




20.2 


± 2.7 


OH 


2 n 3/2 - 2 n 1/2 J= 3/2+ - 5/2- 


96.314 




0.88 


± 





05 




10.5 


± 0.6 


OH 


2 n 3/2 - 2 n 1/2 j= 3/2- - 5/2+ 


96.368 




0.74 


± 





03 




13.4 


± 0.6 


CO 


27 - 26 


96.773 




3.82 


± 





11 




0.2 


± 0.3 


OH 


2 °l/2 J = 5/2~ - 3/2+ 


98.725 




3.61 


± 





33 




38.8 


± 3.5 


OH< 2 > 


2 n 1/2 J= 5/2+ - 3/2- 


98.737 














15.7 


± 3.5 


o-H 2 


505 - 4 14 


99.492 




5.61 


± 





12 




22.2 


± 0.5 


H^O' 1 ' 


5o5 — 4l4 


99.784 




0.03 


±0.01 




24.6 


± 0.8 


CO 


26 - 25 


100.461 




3.44 


± 





00 




11.5 


± 0.2 


P-H3O+W 


27 - if 


100.577 




0.03 


± 





01 








o-H 3 0+' 1 ^ 


2- - 1 + 


100.869 




0.09 


± 





05 








o-H 2 


5i4 — 


100.913 




3.05 


± 





17 




21.1 


± 1.2 


P-H 2 




100.983 




2.35 


± 





00 




29.9 


± 1.2 


p-H 2 


— ^33 


101.210 




0.19 


± 





04 




14.8 


± 2.9 


H 2 S 




102.008 


0.22 ± 0.03 










-11.7 ± 1.3 






P-H 2 


642 - 633 


103.916 




0.08 


± 





02 




19.4 


± 5.6 


P-H 2 


615 - 606 


103.939 




0.11 


± 





02 




18.5 


± 3.7 


o-H 2 


634 - 6 25 


104.090 




0.61 


± 





11 




24.2 


± 4.3 


CO 


25 - 24 


104.445 




3.67 


± 





10 




8.32 


± 0.2 


o-H 2 


2 2 i - lio 


108.073 




3.22 


± 





08 




29.1 


± 0.7 


CO 


24 - 23 


108.763 




3.28 


± 





12 




9.7 ± 0.4 


H^O 


2 2 l - lio 


109.350 


0.08 ± 0.02 


0.18 ± 


0.013 


-8.1 ± 1.7 


48.9 ± 3.5 


p-H 2 


524 - 5i5 


111.626 




0.43 


± 





03 




26.7 ± 1.7 




7 43 - 7 34 






0.09 


± 





02 




11.6 ± 2.7 


CO 


23 — 22 


113.458 




4.82 


± 





15 




10.4 ± 6.3 


o-H 2 


4l4 — 3o3 


113.537 




5.10 


± 





13 




22.6 ± 0.6 


p-H 2 


533 5 24 


113.944 




1.38 


± 





06 




27.6 ±1.3 




^14 3q3 


114.297 




0.16 


± 





01 




44.6 ± 1.7 


OH 


2 n^ ,r, - 2 TT 1 T— 7/2— - 5/2 + 
3/2 L1 l/2 J — '/-* J / * 


115.153 




0.18 


± 





01 




20.2 ± 1.6 


CO 


22 — 21 


118 581 




6.58 


± 





00 




12.5 ± 0.2 


OH 


2n , T — 5/9— '1/9 + 
ii 3/2 J ~ D / z a l z 


119 234 


n °.q -4- n n7 
u.oy in u .u / 


1.70 


± 





17 


-45 5 ± 8 4 


33.9 ± 3.3 


OH 


il 3/2 J — 5/2^ - 3/2 


119.442 




2.42 


± 





16 




35.5 ± 2.2 


18 OH< 1 ' 


2 n 3/2 j= 5/2+ - 3/2- 


119.966 


0.11 ± 0.04 


0.10 ± 0.025 


-25.8 ± 9.1 


51.4 


± 12.8 


18 OH (l) 


2 °3/2 J = 5/2~ - 3/2+ 


120.172 


0.03 ± 0.01 


0.031 


±0.005 


-5.1 ± 0.9 


28.9 ± 4.6 


o-H 2 


432 - 4 2 3 


121.719 




2.28 


± 0.13 




28.1 


± 1.6 


[Nil]' 1 ) 


3 Pl - 3 P 2 


121.890 




0.93 ± 0.47 








CO 


21 - 20 


124.193 




7.89 


± 0.58 




15.0 ± 1.1 


0-NH3 


4 3 - 3 t 


124.647 




0.08 


± 


0.01 




40.8 


± 2.3 


p-NH 3 (1) 


4 , - 3+ 
2 2 


124.796 




0.03 


± 


0.01 




11.5 ± 0.5 


0-NH3 


4 - 3+ 


124.913 




0.05 


± 


0.01 




8,0 


± 0.6 


HDO 




124.954 




0.02 ± 0.004 




9.1 


± 1.9 


p-H 2 


4^4 - 


125.354 




2.25 


± 





00 




25.5 


± 1.1 


P-H 2 


3-J1 - 322 


126.713 




0.62 


± 





08 




30.5 


± 4.1 


o-H 2 


7 25 - 7 16 


127.882 




0.24 


± 





03 




26.4 


± 3.1 




21 - 20 


129.891 




0.02 


± 





01 




10.3 


± 3.9 


CO 


20 - 19 


130.369 




6.19 


± 





22 




16.1 


± 0.6 


o-H 2 


5l4 - 5o5 


134.935 




0.86 


± 





05 




25.4 


± 1.6 


HDO 


4 3i - 4 22 


135.425 




0.08 


± 





02 




17.5 


± 4.1 


o-H 2 


3 30 - 3 2 i 


136.494 




0.71 


± 





05 




30.9 


± 2.3 


CO 


19 - 18 


137.196 




5.22 


± 





05 




18.4 


± 0.2 


P-H 2 


3i3 - 2 02 


138.527 




2.67 


± 





03 




26.4 


± 0.3 


P-H 2 


4i3 - 3 22 


144.517 




1.15 


± 





07 




18.4 


± 0.2 


CO 


18 - 17 


144.784 




3.98 


± 





26 




15.5 


± 1.1 


[OI] 


3 P - 3 P 1 


145.525 




1.17 


± 





19 




28.5 


± 4.6 


P-H 2 


431 - 4 22 


146.919 




1.14 


± 





00 




22.1 


± 1.1 


CO 


17 - 16 


153.267 




8.66 


± 





13 




14.6 


± 0.2 


P-H 2 


3 22 - 3i3 


156.193 




1.23 


± 





08 




32.5 


± 2.2 


o-H 2 


5 2 3 - 4 3 2 


156.266 




1.09 


± 





22 




16.7 


± 3.4 


[CII] 


2 P 3 /2 - 2 Pl/2 


157.741 




2.23 


± 





06 




14.9 


± 0.4 


P-H 2 


532 - 5 2 3 


160.504 




0.45 


± 





00 




24.9 


± 3.2 


CO 


16 - 15 


162.812 




9.48 


± 





20 




12.7 


± 0.4 



Tentative detections Blended lines L01 detection 



10 M.R. Lerate et al. 



Table 4 - Continued 



Species 


Transition 


Rest 
Wavelength 

((im) 


Flux 
absorption 

(10" 17 W cm" 2 ) 


Flux 

emission 
(10- 17 W cm- 2 ) 


Peak Peak 
absorption emission 
(km s -1 ) (km s _1 ) 


OH 


n 3/2 J=3/2+ - 1/2 


163.131 




0.81 ± 0.19 


17.4 ± 4.1 


OH 


2 n 3/2 .1=3/2- - 1/2+ 


163.397 




1.66 ± 0.38 


17.7 ± 4.1 


p-NH 3 


3 2 - 2 t 


165.596 




0.11 ± 0.02 


14.4 ± 2.3 


p-h 2 


633 - 624 


170.131 




0.27 ± 0.01 


43.5 ± 1.3 


co< 3 > 


15 - 14 


173.631 




9.80 ± 0.69 




H 2 8 ° 


3()3 _ 2 12 


174.379 




1.01 ± 0.06 


3.2 ± 1.1 


H^of 1 ' 


303 _ 2 12 


174.491 




0.01 ± 0.004 


4.8 ± 1.9 


o-H 2 


303 _ 2 12 


174.626 




2.38 ± 0.19 


20.1 ± 1.6 


o-H 2 


4 32 - 5o5 


174.914 




0.16 ± 0.01 


23.5 ± 1.3 


o-H 2 


2 12 - loi 


179.527 




2.55 ± 0.31 


23.8 ± 2.9 


o-H 2 


2 21 - 2 12 


180.488 




1.03 ± 0.11 


23.8 ± 2.5 


p-H^O+f 2 ' 


h - 4 


181.054 




0.05 ± 0.01 


11.6 ± 0.3 


h| 8 o< 2 ' 


2 21 - loi 


181.053 








c q(3)(2,3) 


14 - 13 


185.999 




7.58 ± 0.22 




p-H 2 0< 2 - 3 > 


4l3 - 4 4 


187.110 









Tentative detections Blended lines ^ LOI detection 
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Figure 6. The complete Orion KL spectral survey obtained by the ISO LWS in Fabry-Perot modes L03 and L04. The grating observation 
(L01 mode) is superimposed in regions where high resolution data are not available. 
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Figure 6. Continued. Note that features with no label on the L01 spectrum are non real features transferred by the calibration stage. 
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Figure 7. Transitions of [O in], [N III], [O i] and [C II] detected 
by the ISO LWS spectral survey towards Orion KL. Note that 
the LSR velocity of the line peaks (except for the [O i] line at 
145.5 /im) is 10 km s -1 , in agreement with the velocity of the 
quiescent gas (si 9 km s _1 ). 
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Figure 8. Tentative detection of the 122 fim line of [N II]. The 
broadness is due to hyperfine components in the transition as 
can be seen in the detection of the [N n] 122 (im line in Sgr B2 
(Polehampton et al. 2006 in prep.) 



These lines were detected by the Short Wavelength 
Spectrometer (SWS) towards IRc2 and the inferred [S ill] 
33.5/18.7 fim ratio was found to be ~ 0.55 (van Dishoeck 
et al. 1998). This ratio gives an 7V e (S III) of « 1950 cm -3 , 
which exceeds N e (0 ill) by a factor of two, probably due to 
the higher critical densities of the [S ill] lines. 



4.2 [O i] and [C n] 

The [C II] 157.7 ^m and [O i] 63.2, 145.5 fim lines are 
important coolants in photodissociation regions (PDRs), 
whose heating is thought to be dominated by energetic 
photoelectrons ejected from dust grains following FUV pho- 
ton absorption (Tielens & Hollenbach 1985). As important 
coolants, these lines can also originate in shocks driven by 
jets and outflows from young stellar objects (Hollenbach & 
McKee 1989). 



The Orion nebula is an H II region which is ionised by 
the Trapezium stars. The luminosity of the ionising stars 
is about 1O J L0, located approximately 0.15pc from the 
molecular cloud leading to an incident FUV flux of about G 
« 10 5 G (Tielens & Hollenbach 1985). The fraction of the 
FUV flux converted into line emission in these transitions 
is generally of the order 10~ 2 - 10~ 3 (Tielens & Hollenbach 
1985) and is a function of both the gas temperature and 
the electron density. 

The [C II] 157.7 /[O i] 63.2 ^m ratio in Orion KL in- 
ferred from our observations is 3.42 x 10" , in good 
agreement with Tielens & Hollenbach's predictions. The 
[O i] 63.2/145.5 /an ratio combined with the [O i] 63.2 
fim/[C n] 157.7 fira flux ratio can also be used to derive 
the temperature and density of PDRs (Watson 1983). 
Using the PDR temperature density plane defined in Liu 
et al. 2001, we inferred a temperature of 300 K and a 
density of log Nh=5A (cm" 3 ). Considering this, it seems 
natural to postulate that the lines originate in the PDR 
region. However, Orion KL could also be an important 
shock region (Chernoff, Hollenbach & McKee 1982) and 
a potential major contribution to [O i] emission from 
dissociated shocked gas should be taken into account. 
Table |S] lists the surface brightnesses measured with the 
80" LWS beam compared with those observed towards 1 
Ori C (column 3) by the KAO with a « 45" beam (Melnick, 
Gull, & Harwit 1979; Ellis & Werner 1985; Werner et al. 
1984; Stacey, Smyers, Kurtz & Harwit 1983; Russell et al. 
1980; Phillips & Huggins 1981). Column 5 lists predicted 
line surface brightnesses from the photodissociation model 
of Tielens et al. (1985) and column 6 lists the predicted line 
surface brightnesses from the dissociative shock model of 
Hollenbach & McKee (1989) (for n = 10 5 cm" 3 and v = 
30 - 80 km s" 1 ). 

We found that although our [O i] 63.18 fim and [O i] 145.52 
/im line intensities are in good agreement with the shock 
model prediction, the observed [C n] line intensity is much 
higher than predicted by shock models. A PDR model can 
reproduce the ISO LWS [O i] 63.2 fira and [C n] surface 
brightness levels within 35%, although over-predicting the 
LWS [O i] 145.5 lira emission by a factor of 2.7. We note, 
however, that for the smaller beam KAO observations the 
PDR model can match the observed surface brightness 
levels for all three lines. 



4.3 Molecular species 

4-3.1 Water lines 

High velocity gas was first detected at the centre of the 
Orion-KL region as broad wings on 'thermal' molecular lines 
in the millimeter range and as high velocity maser features 
in the 22 GHz line of H 2 (Genzel et al. 1981). These high 
velocity motions may be caused by mass outflows from newly 
formed stars. Many theoretical studies of the Orion region 
(e.g. Draine & Roberge 1982; Chernoff et all. 1982; Neufeld 
& Melnick 1987) have concluded that the rich emission spec- 
trum from thermally excited water vapour should play a sig- 
nificant role in cooling the gas. 

In conditions where the temperature exceeds ~ 400 K, OH 
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Table 5. Detected fine-structure lines 



Species 


Transition 


Rest Wavelength (/^m) 


Flux (W cm" 2 ) 


[O in] 


3 Po - 3 Pl 


51.814 


(3.71 ± 0.005)xl0" 16 


[N in] 


2 P3/2 ~ 2 Pl/2 


57.329 


(5.97 ± 0.12)xl0~ 17 


[O i] 


3 Pl - 3 P 2 


63.184 


(6.53 ± 0.05)xlCT 16 


[O in] 


3 Pi- 3 Po 


88.356 


(1.52 ± 0.03)xlO- 1G 


[Nil] 


3 Pl - 3 P 2 


121.889 


(9.3 ± 4.7)xl0" 19 


[O i] 


3 Pi - 3 Po 


145.525 


(1.17 ± 0.19)xl(r 18 


[C ii] 


2 P3/2 ~ 2 Pl/2 


157.741 


(2.23 ± 0.06)xlCT 19 



Species 


Wavelength (/im) 


KAO at 6 1 Ori C 
(ergs cm -2 s" 1 sr" 1 ) 


ISO LWS BN/KL 

(ergs cm -2 s _1 sr~ 


PDR model 
(ergs cm -2 s _1 sr _1 ) 


Shock model 
(ergs cm -2 s~ 1 sr _1 ) 


[Oi] 


63.18 


4-6xlCT 2 


5.5X10" 2 


4.6xl0~ 2 


6.0X10" 2 


[Oi] 


145.52 


3-6 xl(T 3 


1.3xl0~ 3 


3.5X10" 3 


1.2xl0 -3 


[C ii] 


157.74 


4-7xl0 -3 


2.5xl0" 3 


3.8xl0" 3 


< lxKT 4 



Table 6. Comparison of the observed line intensities towards Orion KL with the PDR model calculations of Tielens et al. (1984) and 
the shock model of Hollenbach et al. (1989). 



and H2O are rapidly formed through the following reactions 
(Elitzur & Watson 1978): 

+ H 2 ^OH + H (1) 

OH + H 2 ^ H 2 + H (2) 

At higher temperatures, OH is also formed via molecu- 
lar oxygen destruction: 

2 + H 2 ^ OH + OH (3) 

eventually being processed into water again by reaction (2). 
Consequently, OH and H2O fractional abundances increase 
in shocks and so does the overall rate of reaction. 
Molecules and atoms behind a shock front act as cooling 
agents by emitting infrared and microwave photons. Colli- 
sions between molecules at high temperature populate the 
vibrational and rotational levels of the molecules. Following 
these reactions, rotational and vibrational water transitions 
are predicted to occur, providing a ubiquitous tracer of 
shock-heated gas. 

The widespread nature of the water vapour around IRc2 
has been probed with maps at 183 GHz (Cernicharo et al. 
1990, 1994, Cernicharo & Crovisier, 2005), the first time 
that its abundance was estimated in the different large-scale 
components of Orion IRc2. Harwit et al. (1998) analysed 8 
lines of water observed by LWS FP in L04 mode, concluding 
that these lines arise from a molecular cloud subjected to a 
magnetohydrodynamic C-type shock. From their modelling, 
they derived an H2O/H abundance of 5 x 10 -4 . However, 
the interpretation of these lines in the ~ 80" LWS beam and 
the determination of the water abundance in the different 
components of Orion remains a long standing problem. This 
is basically due to two main issues: the complexity of the 
different dynamical and chemical processes that take place 
within the region encompassed by the LWS beam, including 
outflows and several gas components, and the need for new 
H2O collisional rates appropriate for the temperatures pre- 
vailing in shocks. 

Water lines appear in the survey as resolved (typically 70 



km s _1 FWHM) with a total of 70 detected lines. The 
line profiles range from predominantly P-Cygni at shorter 
wavelengths to predominantly pure emission at longer wave- 
lengths (see Figure [§J . Radial velocities appear to be cen- 
tred at ~ - 15 km s _1 in all absorption lines (shorter wave- 
lengths) consistent with the results found in the SWS range 
(from 2 fim to 45 ^m, van Dishoeck et al. 1998; Wright et 
al. 2000). However, the pure emission lines of H2O peak at 
around ~ +30 km s _1 (see Section 4.4), where the velocity 
of the quiescent gas is 9 km s _1 (Cohen et al. 2006). The 
same behaviour is found for OH lines, where this has been 
interpreted as evidence of an outflow at a velocity of 25 
± 5 km s _1 (see discussion in Section 4.3.3; also Goicoechea 
et al. 2006). 

4-3.2 Water isotopes 

Several transitions of the water isotopic variants H^O, 
H^O and HDO are detected in the survey (Figure IT01 and 
Figure ITTI . The detection of isotopes is of special interest 
for abundance determinations and the interpretation of the 
spatial origin of the lines. These lines are excellent water 
tracers as they are more likely to be optically thin, so that 
they can be used to derive the water abundance via the 
adoption of a [ ie O]/[ 18 0] isotopic ratio. 
The detections of H^ 8 at 67.19 /an and 102.01 /im are 
of particular interest (see Figure HOB . These H3, 8 lines 
show a main absorption component that contrasts with 
the equivalent H^O transitions at 66.43 /im and 100.98 
(im, which appear as a P-Cygni and an emission line, 
respectively. The fact that both lines peak at ~ - 10 km 
s _1 LSR can be interpreted as indicating an origin in the 
main outflow. 

The H3> 8 rotational diagram (Figure [T21 leads to a column 
density of iV coi « 2.8 x 10 14 cm' 2 and T rot = 60 K. If we 
adopt [ 16 0]/[ ls O] = 500 (Bergin et al. 1998), the estimated 
total H2O column density is ~ 1.5 x 10 17 cm -2 . 
However, the interpretation of the exact spatial origin of 
the water lines is not straightforward and requires more 
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Figure 9. Water lines detected by the survey towards Orion KL 
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Figure 10. A selection of water isotope lines detected by the 
survey 



Figure 11. Tentative detections of five transitions of HDO ob- 
served towards Orion KL 



sophisticated models. A detailed study of these lines will be 
published in future papers (Cernicharo et al. in preparation) 
and is beyond the scope of this work. 

A total of 5 far-infrared HDO lines are tentatively iden- 
tified for the first time towards Orion KL, from 330 — > 2n 
at 59.93 (im to 4 3 i -> 4 22 at 135.42 /Ltm (see Figure ITU. The 
HDO rotational diagram gives a total column density of N co i 
~ 3 x 10 14 cm~ 2 . This value is significantly lower than that 
obtained from millimeter wave observations (Turner et al. 
1975, Beckman et al. 1982, Moore et al. 1986). However, if 
the main HDO emission originated in the Hot Core, a low 
column density could be explained if the emission from the 
core is blocked by the intermediate-velocity gas. This hy- 
pothesis was previously pointed out by Pardo et al. (2001) 
in their analysis of the first detections of two submillimeter 
HDO lines from the KL region (transitions 2i2 — * In at 
848.9 GHz and In -»• Ooo at 893.6 GHz), where the total 
HDO column density was estimated to be in the range (4 - 
6) x 10 16 cm -2 . 

They concluded that the HDO transitions detected at above 
800 GHz arise mainly from a very compact region (HPW 
15") in the Plateau. Our lower column density result (~ 
10 14 cm -2 ) could therefore be explained if the compact HDO 
emission is diluted within the LWS beam. Special care has 
to be taken in the interpretation of tentative detections as 
the rotational diagram for these HDO lines could underes- 
timate the total column density. More sophisticated models 
are needed to aid the interpretation. 

4.3.3 OH Lines 

Far-infrared line emission from excited OH rotational states 
was first detected toward the em bedde d star-forming region 
in Orion KL by Storey et al. <198lfl using the KAO. In 
these early measurements, the 2 n 3 / 2 J=5/2 — > 3/2 features 
at 119.23 lira and 119.44 /im were seen in emission. The 
fact that the excitation temperature above the ground state 
for J up =5/2 is about 121 K and that the OH was observed 
fs 30" north of KL, led these authors to assume that the 
emission they detected came from the shocked gas region 
surrounding BN/KL. 



Watson et al. l)l985l) reported KAO observations of the 
two lowest lying pure rotational transitions in the 2 n 3 / 2 
electron ic sta te of OH, at 119 fim and 84 yum. Melnick 
et al. il990f) modelled these detections, along with their 
own KAO detections of the 2 Ii\/2 J=3/2 — t- 1/2 line at 
163 [im, the 2 Hi /2 -> 2 n 3/2 J=3/2~ -» 3/2+ line at 53 
fim and the 18 OH 2 n 3/2 J=5/2+ -> 3/2" line at 119.44 
/im. They found that the best fit to the data required the 
presence of at least three components (see discussion below). 

A total of 22 transitions of OH are detected by our 
survey, ranging from the 2 ^-i/2 ~ 2 H 3/ / 2 J=5/2~^5/2 + 
transition at 48.7 /jm to the n 3 / 2 J=3/2~^l/2 + tran- 
sition at 163.4 /im. These transitions range up to upper 
energy levels of ~ 618K. Their line profiles show a similar 
behaviour to those of water, evolving from pure absorption 
or P-Cygni at shorter wavelengths to pure emission at 
longer wavelengths. 

Table □ lists a comparison of the fluxes and surface bright- 
nesses of seven OH transitions measured with different 
beam sizes by different observers. 

A similar behaviour to that of the water line velocities is 
also found for the OH detections, showing center velocities 
between - 10 km s _1 and +30 km s _1 (see Section 4.4) 
when the line appear in absorption or emission respectively. 
Radiative transfer modelling has been recently performed 
for the OH lines (Goicoechea et al. 2006) and concludes that 
most of the newly detected excited OH lines in the survey 
originate in a compact region (D « 25"), not resolved by the 
LWS, that forms part of the low velocity plateau component. 
Average values of the physical conditions in this region are 
estimated to be: n(H 2 ) w 5 x 10 5 cm~ 3 , T k w 100 K and 
an abundance of x(OH) w (0.5 - 1) X 10" 6 . 

4-3.4 OH isotopes 

The detection of 18 OH at the sensitivity of the FP-LWS ob- 
servations was found to be very difficult. The weakness of 
the lines makes them easily confused with noise and only 
the 18 OH 2 n 3/2 5/2+^3/2~ transition at 119.96 /im and 
the 2 n 3 / 2 5/2~— >3/2 + transition at 120.17 /xm were tenta- 
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Figure 12. Rotational diagrams i.e. the natural logarithm of the column density in the ith stage, Ni, divided by the degeneracy gi versus 
the upper energy E up for the HDO and H 2 s O transitions detected by the survey towards Orion KL. Due to the high uncertainties of the 
high energy HDO lines, only the lowest energy transitions are considered. Note also that P-Cygni lines are excluded from the diagram. 



Line 



Wavelength 
H 



FWHM Beam Size 
(arcsec) 



Reference^ 



(10 



Flux 
' 17 W cm" 



2)(b) 



Surface Br. 
-2 



(10 3 ergs cm 2 s 1 sr 1 ) 



2 n 1/2 j= 3/2- 

2 n 1/2 j= 3/2- 1 

2 n 3/2 J= 5/2^ 
2 n 3/2 j= 5/2- 



2 n 3/2 J= 7/2^ 
2 n 1/2 j= 3/2- 



+ 1/2+ 

♦1/2- 

+3/2- 
+3/2+ 



2 n 3/2 J= 7/2--^5/2+ 



+5/2- 

+ 2 n 3/2 J=3/2+ 



163.396 

163.121 

119.441 
119.234 

84.597 

84.420 
53.351 



67 
55 
60 
67 
55 
60 
78 
Id 
78 
44 
45 
77 

30 
60 
77 
30 
60 
85 

40 



1.66 ± 0.38 
1.3 
1.2 ± 0.6 
0.81 ± 0.19 
1.3 
1.2 ± 0.6 
2.42 ± 0.13 
0.7 ± 0.3 
1.76 ± 0.17 
0.83 
1.88 

(1.13 ± 0.38) ahs 
(1.34 ± 0.26) em 
0.56 
1.4 ± 0.4 
1.72 ± 0.13 
<0.5 
1.0 ± 0.3 
(2.45 ± 0.26) abs 
(0.26 ± 0.13) em 
3.00 abs 
0.63 em 



2.04 

1.6 

1.2 

0.99 

1.6 

1.2 

2.16 

1.3 ± 0.5 

1.57 

1.5 

3.3 

1.01 abs 

1.19 em 

2.1 

1.4 

1.54 

<1.9 

1.0 

1.86 afts 

0.19™ 
66 abs 

1.4 em 



( a ) References : 

1 Values using ISO LWS 
data in this work 

2 Melnick et al. 1987 

3 Viscuso et al. 1985 

4 Watson et al. 1985 

5 Melnick et al. 1990 

W 

(abs) p Cygn j 

absorption component 

(em) p Cygn j 

emission component 



Table 7. Comparison of OH line measurements obtained with different beam sizes and instruments. The uncertainty in the absolute 
fluxes and intensities is ± 30% unless indicated otherwise 
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Figure 13. OH lines observed in the survey 
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tively identified (see Figure [Till . 

As with water, isotopic variants play an important role in 
the interpretation of the OH lines. Both 18 OH transitions de- 
tected show a difference in their line profiles when compared 
with their corresponding ie OH transitions. As the 18 OH op- 
tical depths are smaller than those measured for the 16 OH 
transitions, and since both line profiles are P-Cygni, the 
18 OH lines indicate that OH is predominantly associated 
with the Plateau outflow (see also Goicoechea et al. 2006). 

4-3.5 CO emission lines 

High-J CO transitions are well known among the most 
common tracers of highly excited regions. OH and H2O can 
react with C + via the reactions: 



C + + OH -» CO + + H 
C + + H 2 -> HCO + + H 
CO + + H 2 ^ HCO + + H 
HCO + +e^CO + H 



(4) 
(5) 
(6) 
(<) 



These reactions are sufficiently rapid to ensure that sub- 
stantial amounts of CO are produced in shocks (Pineau des 
Forets et al. 19 86, 19 87) 

Watson et al. I|198fl) made the first far-infrared detection 
of interstellar CO with observations of the J=21 - 20 and 
J=22 - 21 transitions at 124 /jm and 119 /im respectively, 



from the BN/KL region. Using col lision al excitation cross 
sections calculated by Storey et al. (ll98lT) . they derived CO 
rotational level populations for a number of temperatures 
and densities. They showed that the observed line intensities 
could be modelled by emission from two components: a 
2000 K component with n(H 2 )~l x 10 6 cm" 3 and a 400 
- 1000 K component with n(H 2 )~5-2 x 10 6 cm" 3 . More 
recent studies using ISO data indicate that the CO emission 
can be modelled using three temperature components, 
describing the plateau and the ridge emission (Sempere et 
al. 2000, Maret et al. 2001) 

CO detections originating form J up =14 to J up =52 are 
identified by our survey in a total of 26 detected emission 
lines. Line profiles are shown in Figure Upl and the measured 
line fluxes and intensities are com pared in Table |H| with 
th ose m easured by Storey et al. il98lt) and Watson et 
al. Jl985l) using the KAO. Note that the field of view for 
KAO's observations was 60 arcsec for the i up = 17, 16 lines 
and ~44 arcsec for the 3 up — 21, 22, 26, 27, 30 and 31 lines. 
An interesting result from Table |S| is that almost identical 
surface brightnesses were measured for the J=17-16 and 
J=16-15 transitions when observed with the 60" KAO 
beam and the 80" LWS beam. However for the higher J 
CO lines we find that the surface brightness in the 44" 
KAO beam is about 60% higher than in the 80" LWS 
beam, while the integrated fluxes in the two beam sizes are 
similar, indicating that the source size may be smaller than 
44 arcsec for these transitions. 

CO is also an important tracer of H 2 , since it is the second 
most abundant molecule in the interstellar medium after 
molecular hydrogen. Despite theoretical and observational 
uncertainties in the use of a canonical iV(CO)/iV(H 2 ) ratio 
(Williams et al. 1984; van Dishoeck et al. 1992; Sakamoto et 
al. 1996), theoretical studies of the CO/H 2 abundance ratio 
have concluded that on large scales it can be considered to 
be constant (Taylor et al. 1993). 

In Figure [TrU we show rotational diagrams for three dif- 
ferent CO J-ranges, finding three different rotational tem- 
peratures. Belo w, we compare our results with the model of 
Sempere et al. a2000f) which was developed using a radiative 
transfer model fit to ISO data obtained at low spectral res- 
olution (grating mode), together with selected observations 
at high spectral resolution. In this model: 
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Figure 15. CO lines detected by the survey 
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• The CO emission from J=18 to 33 can be explained by 
a two temperature component model of the plateau region; 
the inner region reproduces the emission from J=33 to 28 
and the colder gas contributes to the lower J lines. The high 
J transitions (J> 34) reveal the presence of a very hot gas 
component ( T w 1500-2000K). 

• They assumed densities of 10 7 cm -3 and a temperature 
of 400 K for the inner region and 10 6 cm -3 , 300K for the 
external part of the plateau. 

• The resulting column densities in Sempere et al's model 
were: iV(CO) = 10 19 cm" 2 and iV(CO)=3.5 x 10 18 cm" 2 for 
the inner and outer plateau, and 7V(CO)=10 17 cm -2 for the 
hot gas component. 

The results from our rotational diagrams shows consistency 
with the different temperature components predicted by 
Sempere et al's model. Our results confirm the presence of 
hot gas (T « 660 K; Figure [THJ which differs by w 300 K 
from the warm component (T ~ 360 K, see Fieure lTrjl . How- 
ever, we found that this warm component has a higher CO 
column density than that found by Sempere et al., with a 
column density of the order of ~ 2.3 x 10 19 cm -2 . Accord- 
ing to our results, the inner part of the Plateau component 
is traced by the warm gas which emits the CO transitions 
with J < 28 with a rotational temperature of ~ 360 K. 



4-3.6 CO isotopic variants 

CO isotopic variants are hardly detected in the survey. Only 
one line of 13 CO (the J = 21-20 transition at 129.89 (im) is 
tentatively detected with a flux of 2 ± lxlO -19 W cm -2 . 



4.4 Line kinematics 

Although kinematical properties are diluted in the large 
LWS beam, resolved line profiles and velocity peaks can 
trace the overall dynamical gas properties. Figure H71 shows 
the upper transition energy of the main molecular detec- 
tions of H2O, OH and CO, as a function of the emission and 
absorption LSR line velocity peaks. For the H2O and OH 
lines, both pure emission or absorption and P-Cygni peaks 
are plotted. The velocity of the quiescent gas is 9.0 ± 0.5 
km s" 1 (Cohen et al. 2006). Subtracting this value from the 
velocity peaks gives an indication of the exact blue- or red- 
shift of the lines. Several conclusions can be deduced from 
Figure 1171 

• H2O and OH radial velocities trace the velocity dis- 
tribution of the same expanding gas. Absorption lines are 
centred between -15 and -30 km s _1 and emission lines be- 
tween + 15 and +45 km s - . This confirms previous analy- 
ses of these lines (Harwit et al. 1998; Cernicharo et al. 1999; 
Goicoechea et al. 2006) which interpreted them as evidence 
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Figure 16. Rotational diagrams i.e. the natural logarithm of the column density in the ith stage, Ni, divided by the degeneracy gi 
versus the upper energy E up for the observed CO emission towards Orion KL using ISO LWS-FP data. 
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(a) Line fluxes and intensities averaged over a 80" beam observed by the ISO LWS in Fabry-Perot mode. 

(ft) Line fluxes and intensities averaged over a 44" beam measured by Watson et al. (1981) using the UC Berkeley tandem Fabry-Perot 
spectrometer. 

(c) Storey et al. (1981) using a 60" beam. 
Table 8. LWS Fabry-Perot CO detections towards Orion KL, compared with those of Watson et al. (1981) 



for outflows driven by the star formation activity of Orion 
KL. 

• For H2O and OH, the absorption and emission com- 
ponents of the P-Cygni lines tend to be at more extreme 
velocities than those of the 'pure' absorption or emission 
lines. Detailed radiative transfer modelling is underway to 
attempt to match this behaviour. 

• The CO radial velocities indicate a different excitation 
mechanism to that of H2O and OH, with emission peaks 
centred between (5 - 15) km s _1 . Considering an average 
velocity uncertainty of ±10 km s , this range is consistent 
with the velocity of the quiescent gas (w 9 km s _1 ). As 
very high energy CO transitions are detected, the lines may 
originate in the hot and warm quiescent post-shocked gas. 

4.5 HsO + line detections 

HsO + is one of the key species in the interstellar chemistry 
of oxygen. This saturated molecular ion leads to the forma- 
tion of OH and H2O by dissociative recombination, with an 
imprecisely known branching ratio (Bates 1986; Sternberg 
1995). It can be formed from H3" by the reaction: 

#3+ + H 2 -> H s O + + H 2 (8) 

or via the reactions: 



H 3 + + O -> OH + + H 2 (9) 

OH+ + H 2 ^ H 2 + + H (10) 

H 2 + + H 2 ^ H a O + + H (11) 

Its formation is more likely to occur via reactions 9- 
11 (Wootten et al. 1986). Destruction of HsO + occurs via 
electron recombination: 

H s O + + e~ -> H 2 + H (12) 

H s O + + e~ -> OH + H 2 (13) 

leading to the production of H2O and OH. 
HsO" 1 " has a pyramidal structure; inversion transitions are 
produced when the oxygen atom tunnels through the plane 
of the hydrogen atoms. The ground state inversion splitting 
is « 55 cm -1 (Liu et al. 1985). This large value makes the 
fundamental transitions of the v 2 mode lie at submillimeter 
and far-infrared wavelengths (Wootten et al. 1986, Bogey et 
al. 1985). H 3 0+ was first detected at 365 GHz in OMC-1 and 
Sgr B2 by Wootten et al. (1991) and more recently in the 
far-infrared towards Sgr B2 (Goicoechea and Cernicharo, 
2001; Polehampton et al. 2006). Wootten et al. (1991) 
modeled the excitation of the para-HsO + transition at 365 
GHz, finding abundances x(H 3 0+) six 10~ 9 -> 5 x 10" 9 . 
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Figure 17. LSR peak velocities of H2O, OH and CO lines observed towards Orion KL with the ISO LWS-FPs. Note that the LSR 
velocity of the quiescent gas is fa 9 km s _1 LSR (Cohen et al. 2006). Errors are from goodness of fit estimations. 



The lowest rotational levels of the ortho and para lad- 
ders of the V2 ground-state H30 + inversion mode (0 + — >0 _ ) 
are detected in the survey. These are the rotational inver- 
sions 2^^1^ at 100.58 /Ltm and 2^— >lj at 100.87 /im and 
the pure inversion lj"— >1^" transition at 181.05 fim (see Fig- 



ure I18t . The optical depths of these transitions estimated 
with the RADEX code are r > 10. We used as input the 
column density values inferred by Wootten et al. (1991; fa 
10 14 cm" 2 ), T K =80 K and ?ih 2 = 10 5 cm" 3 . The rotational 
diagram method cannot therefore be applied to derive phys- 
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Figure 18. H30 + lines observed towards Orion KL 

ical parameters. Further modelling of optically thick lines 
including radiative excitation by the FIR dust continuum 
emission is needed and will be published in forthcoming pa- 
pers. 

4.6 NH3 line detections 

Like H30 + , the NH3 molecule has a pyramidal structure. 
Inversion transitions are produced when the nitrogen atom 
tunnels through the plane of the hydrogen atoms (analogous 
to the tunnelling of the O atom in H3O 4 "), with the signifi- 
cant difference that the inversion splitting m 1 cm -1 , much 
smaller than that of H 3 0+ (Ho & Townes, 1983). 
NH3 was first detected at far-IR wavelengths by the KAO 
towards Orion-KL (Townes et al. 1983), in the rotational 
1/2=1; - 3^~ line at 124.6 fim. The fact that this line is 
very optically thick (r « 10 3 ) and that it is seen in emission, 
led the authors to consider that this level is collisionally ex- 
cited, probably originating in the Hot Core component. 
The 4J - 3jj~ transition at 124.6 /xm is also detected by the 
LWS, as well as the following transitions: 4JT - 3^ 124.8 /im, 



4q - 3+ 124.9 fim and 3? - 2+ 165.6 /im (see Figure HH. 



We estimated the optical depths of these transitions with the 
RADEX code, using as input the inferred column density of 
10 14 cm~ 2 , T X =80K and n H2 =10 5 cm~ 3 . For these condi- 
tions, only the 3JT - 2^~ transition at 165.6 (im is moderately 
optically thick (r = 1.4). Figure |2"U1 shows the rotational di- 
agram of NH3, with an inferred column density of 1.3 x 
10 14 cm -2 and a rotational temperature of 40 K. However, 
if the NH3 column density is ^ 10 15 cm -2 the code predicts 
high optical depths and our rotational calculation could be 
underestimated. Estimations of the NH3 abundance using 
more sophisticated models will be published in forthcoming 
papers. 



4.7 Unidentified and weak features 

Some isotopic lines that were detected by the LWS in L04 
mode were very difficult to detect in the L03 scans, due to 
differences between repeated scans. This affected the detec- 
tion of weak transitions. However, some unidentified lines 
were detected after carrying out an exhaustive analysis, in 
which the lines were seen in at least two different obser- 
vations plotted together in order to discriminate between 
glitches and spurious features. The observations were also 
analyzed before performing smoothing corrections, to ensure 
that no features were missed in this step. Table [TOl gives a 
list of the unidentified features found in the survey. The line 
fluxes were derived using similar fits to those used for iden- 
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Figure 19. NH3 lines observed towards Orion KL 
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Figure 20. Rotational diagram for the NH3 lines detected to- 
wards Orion KL. 



fined lines. 

Amongst the candidate detections are several rotational 
transitions of the H2O vibrational bending mode v2 (44.09 
and 84.36 /im). There are also many rotational transitions in 
the far-IR of slightly asymmetrical species that could be con- 
tributing to the spectrum but may not have been detected 
at the FP resolution, such as HNO, HNCO or HOCO+ and 
the low-energy bending modes of carbon chains. 



5 CONCLUSIONS 

A high spectral resolution 44-188 (im line survey towards 
Orion KL has been carried out with the ISO LWS in Fabry- 
Perot mode. A total of 152 lines are identified and a fur- 
ther 34 lines remain to be identified. A basic analysis of 
the molecular detections was carried out by deriving rota- 
tional temperatures and column densities (see summary in 
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Table 9. Rotational temperatures and beam averaged column densities N co ;, estimated using LTE rotational diagram method or using 
non-Local radiative transfer models, for the different molecular species identified in the far-infrared survey towards Orion KL with ISO 
LWS-FP data. Column 2 list the method used to estimate the column density; [1] rotational diagram method, [2] non-Local radiative 
transfer models, [3] via the isotopic ratio.' 11 ) From Goicoechea et al. 2006 



Wavelength observed Flux absorption emission 

W (10- 18 W cm" 2 ) 



44.09 ± 7.5 x 10~ 3 1.54 ± 0.98 X 

46.29 ± 5.5 x 10" 4 1.46 ± 0.31 X 

48.48 ± 5.5 x 10~ 4 6.36 ± 1.82 X 

48.75 ± 1.1 x IO- 3 2.58 ± 0.83 X 

52.29 ± 5.7 x 10" 4 4.14± 1.61 X 

52.57 ± 9.5 x 10" 4 6.07 ± 2.01 X 

58.02 ± 8.5 x 10~ 3 1.44 ± 0.38 X 

58.24 ± 5.5 x 10" 3 1.67 ± 0.53 X 

63.78 ± 3.6 x 10~ 4 5.06 ± 1.06 X 

64.34 ± 9.5 x 10- 3 0.87 ± 0.27 X 

66.63 ± 8.5 x 10~ 4 2.03 ± 0.74 X 

67.14 ± 4.5 x 10" 3 1.55 ± 0.75 X 

68.61 ± 8.8 x 10~ 3 1.98 ± 0.47 X 

84.36 ± 5.5 x 10~ 3 1.57 ± 0.93 X 

91.13 ± 7.5 x 10- 3 0.36 ± 0.15 X 

91.97 ± 3.2 x 10" 4 4.08 ± 0.10 X 

101.78 ± 4.1 x 10~ 4 0.31 ± 0.10 X 

104.39 ± 5.6 x 10~ 3 0.68 ± 0.37 X 

108.30 ± 2.2 x 10" 4 0.32 ± 0.05 X 

112.82 ± 5.4 x 10" 3 0.29 ± 0.14 X 

114.96 ± 7.5 x 10~ 3 0.27 ± 0.093 X 

115.68 ± 2.5 x 10~ 3 0.41 ± 0.22 X 

116.27 ± 6.9 x 10~ 4 0.89 ± 0.24 X 

116.60 ± 9.5 x 10~ 3 0.49 ± 0.35 X 
122.37 ± 4.4 x 10~ 4 1.21 ± 0.027 X 

123.80 ± 2.7 x 10~ 4 7.40 ± 0.088 X 

126.64 ± 8.5 x 10- 3 0.19 ± 0.056 X 

131.44 ± 9.5 x 10" 4 1.23 ± 0.45 X 

133.34 ± 3.2 x 10" 3 0.76 ± 0.36 X 

136.61 ± 2.3 x 10" 3 0.69 ± 0.20 X 
137.50 ± 8.4 x 10~ 3 0.99 ± 0.48 X 
139.54 ± 4.8 x 10- 3 1.66 ± 0.065 X 
140.60 ± 8.6 x 10~ 3 0.33 ± 0.15 X 
144.47 ± 5.9 x 10" 4 4.11 ± 0.86 X 



Table 10. List of unidentified and weak features found on the survey towards Orion KL with ISO high resolution spectra Fabry- 
Perot mode. The unidentified line selection criteria was based on the features that were seen in at least two different observations 
plotted together in order to discriminate between glitches and spurious features. The observations were also analysed before performing 
smoothing corrections, to assure that no features were missed in this step. 
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Table |5J and by comparison with previous measurements 
and published models. We found that: 

• The spectrum is dominated by the main molecular 
coolants: H2O, CO and OH together with the forbidden lines 
[O I], [O ill], [N ill] and [C II]. 

• The analysis of the [O 1] and [C 11] fine structure lines 
indicates that a PDR model can reproduce the observed 
[O 1] 63.2 fim and [C 11] 157.7 fim surface brightness levels 
although it over-predicts the [O 1] 145.5 /im emission by a 
factor of 2.7. 

• The water and OH P-Cygni profiles, along with the 
kinematical line properties, confirm that most of the de- 
tected emission is associated with gas expanding in the out- 
flow from the KL cluster. 

• The observed molecular emission is consistent with pre- 
vious models where several main components were clearly 
distinguished; the CO detections confirm the differentiation 
into physically distinct components, with column densities 
ranging from ~ 10 19 to « 10 17 cm -2 and temperatures from 
350 to 650 K. 

• HDO and HsO + are tentatively detected for the first 
time in the far-infrared range towards Orion KL. The de- 
rived HDO column density (« 10 14 cm -2 ) is lower than the 
values (10 16 cm -2 ) obtained at millimeter and submillime- 
ter wavelengths. Due to the high optical depths of the lines, 
more sophisticated non-local radiative transfer models are 
needed to estimate the H3O 4 " and NH3 column densities. 

The 80 arcsec ISO LWS beam size encompasses a large 
range of physical conditions, ranging from quiescent cool 
gas to outflows that alter the chemistry of the region. Con- 
sequently, the exact interpretation of the line survey requires 
modelling. We plan to model the chemistry and dynamics 
of the main components using a coupled chemical-radiative 
transfer code. 
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